THE DERIVE - NEWSLETTER #68 

ISSN 1990-7079 


THE BULLETIN OF THE 



Contents: 


1 Letter of the Editor 

2 Editorial - Preview 

3 User Forum 
Bernhard Kutzler 

5 Optimization Problem with TI-Nspire CAS 

Pat and Carl Leinbach 
1 2 Directing Our Suspicions (AHP) 

Alessandro Perotti 

23 Groebner6 . dfw - Grobner Bases in DERIVE 6 
Josef Bohm 

36 Surfaces from the Newspaper (3) 

37 User Forum(2) 

42 Patterns from India 


December 2007 









D-N-L#68 


Book Shelf - TIME 2008 
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[1] Roman Hasek, Uziti Derive ve vyuce Matematiky, Ceske Budejovice 2007, 

ISBN 978-80-239-9054-6 

[2] Pat und Carl Leinbach/ Josef Bohm, Forensische Mathematik fur den Unterricht, 
bk teachware SR-57, ISBN 978-3-901769-8 

[3] Pavel Pech, Selected Topics of Geometry with Classical vs. Computer Proving, 

World Scientific, Singapore, ISBN 978-981-270-942-4, 

If you know about CAS-related publications then please inform me. Books dealing 
with the use of technology in math education in general are also interesting. Josef 


TIME 2008 Call for Papers 

• 30 April 2008: Submission of abstracts on or before this date 

• 31 May 2008: Acceptance of abstract notification by this date 

• 30 June 2008: Early bird registration on or before this date 

• 1 September 2008: Final date for registration and inclusion of abstract in the programme 

• 10 September 2008: Final programme published on the web-page 

• 22 to 26 September 2008: Conference 

• 31 October 2008: Acceptance notification for selected papers for the SAJS. 

Abstracts can be submitted as e-mail attachments to Josef Bohm ( nojo.boehm@pgv.at ) or 
Michel Beaudin ( Michel.Beaudin@etsmtl.ca ) for the Derive/Nspire part of the conference 
and for the ACDCA part to Bernhard Kutzler ( b@kutzler.com ) or Vlasta Kokol-Voljc 
(vlasta.kokol@uni-mb.si ). 

An online submission method is being developed. E-mail attachments need not be re- 
submitted when the online system is up and running. 

The deadline for submission is 30 April 2008 with acceptance of notification by 3 1 May 
2008. 

You can find more important information about TIME 2008 at 

http : / / time . tut .ac.za 


Many thanks to Steve Joubert and his team in South Africa. 


Download all DNL-Derive- and Tl-files from 

http : //www . austromath . at/dug/ 
http : //www . bk- teachware . com/main . asp?session=37 505 9 
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Dear DUG Members, 

At the end of a year it is always time to give a re- 
port of our activities. You can find a summary of 
DUG-year #17 in the box below. 

Despite the fact that distributing DERIVE was 
ended we can report a lot of activities. 

First of all Fd like to draw again your attention to 
the forthcoming Conference in South Africa. In- 
formation can be found on the website and on our 
Information page. 

This DNL contains less contributions as usual. But 
you will hopefully see that they are very interest- 
ing. 

There is the announced article provided by Bern- 
hard Kutzler who demonstrates the interplay of 
various applications of the TI-Nspire CAS. We 
hope that this contribution will inspire and incite 
other Nspire users to share with us their - good or 
bad - experiences with this CAS. It has now been a 
considerable time on the market and we can expect 
one or the other reaction. 

The second article is a result of a cooperation be- 
tween Pat and Carl Leinbach and me. Pat and Carl 
gave a workshop on Forensic Mathematics at DES- 
TIME 2006 in Dresden. I was so fascinated by their 
examples that I suggested to produce a German 
translation adapted for the use in school, and bk 
teachware promised to publish the booklet. Pat and 
Carl provided additional activities for this collec- 
tion. One of them deals with AHP 


(Analytical Hierachy Process) which is a great 
application of matrix calculation. We, Pat, Carl and 
I, are very grateful for bk-teachware's permission 
to present this chapter - in English - in this DNL. 
At the other hand this might be a mouth watering 
for the German readers to learn more about Foren- 
sic Mathematics (see our Book Shelf). 

In the last DNL I announced an article about Grob- 
ner Bases. Alessandro Perotti sent his paper just in 
time. I added two informative examples, which 
were accepted by Alessandro. Many thanks for 
your reliable cooperation. 

The User Forum shows among other questions a 
fine programming competition. 

Just recently my wife Noor and I spent three weeks 
in India (Rajastan). Many palaces and temples were 
decorated with impressive patterns. I could add 
many pictures to my collection of backgrounds for 
DERIVE. You can find some of them on the last 
page. If you are interested in one or the other pic- 
ture, then let me know. I'll send it. 

Finally we will wish you all a Merry Christmas and 
a Happy and Healthy New Year 2008. We hope to 
meet many of you 2008 in South Africa. 



DUG Report 2007-12-21 

2007 was a sad year for DERIVE enthusiasts because DERIVE was ended in many countries. At the 
same time the DERIVE and CAS-TI User Group was asked by TI for cooperation to make the 
TI-Nspire a successor of DERIVE. DUG members participated at the first “Transition Conference” 
which gave the opportunity to meet the Nspire-team to deposit the wishes of the DERIVIANS for 
future releases of the Nspire with a special impact on didactical issues. 

But there are a couple of good news for 2007, too: 

• We could welcome 37 (! !) new members in the DUG 

• We support organizing TIME 2008 in South Africa 

• We could set up a lot of contacts between DUG-members who had special requests 

• We received high quality contributions for the DNL 

• DUG members were involved in national and international projects 

• DUG members attended national and international conferences giving lectures and workshops 

• 4 revised and 4 new DERIVE & CAS-TI Newsletter were published 

• DUG members reached high positions in universities 
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The DERIVE-NEWSLETTER is the Bulletin 
of the DERIVE & CAS- 77 User Group. It is 
published at least four times a year with a 
contents of 40 pages minimum. The goals 
of the DNL are to enable the exchange of 
experiences made with DERIVE, 77-CAS 
and other CAS as well to create a group to 
discuss the possibilities of new methodical 
and didactical manners in teaching mathe- 
matics. 


Editor: Mag. Josef Bohm 
D'Lust 1, A-3042 Wurmla 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English arti- 
cles will be warmly welcomed nonetheless. 
Your contributions will be edited but not 
assessed. By submitting articles the author 
gives his consent for reprinting it in the 
DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 

Next issue: March 2008 

Deadline 15 February 2008 


Preview: Contributions waiting to be published 

Some simulations of Random Experiments, J. Bohm, AUT, Lorenz Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

ANOVA with DERIVE & Tl, M. R. Phillips, USA 

Financial Mathematics 4, M. R. Phillips 

Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE, J. Bohm 
Henon & Co, J. Bohm 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 

An Interesting Problem with a Triangle, Steiner Point, P. Luke-Rosendahl, GER 

Overcoming Branch & Bound by Simulation, J. Bohm, AUT 

Diophantine Polynomials, D. E. McDougall, Canada 

Graphics World, Currency Change, P. Charland, CAN 

Precise Recurring Decimal Notation, P. Schofield, UK 

Cubics, Quartics - interesting features, T. Koller & J. Bohm 

Logos of Companies as an Inspiration for Math Teaching 

Exciting Surfaces in the FAZ / Pierre Charland 's Graphics Gallery 

BooleanPlots.mth, P. Schofield, UK 

What is hiding in Dr. Pest? B. Grabinger, GER 

Truth Tables on the Tl, M. R. Phillips 

Advanced Regression Routines for the TIs, M. R. Phillips 

Where oh Where is IT? (GPS with CAS), C. & P. Leinbach, USA 

Embroidery Patterns, H. Ludwig, GER 

Inspheres and Exspheres of a Tetrahedron, D. Sjostrand, SWE 

Mandelbrot and Newton with DERIVE, Roman Hasek, CZ 

Snail-shells, PiotrTrebisz, GER 

A Conics-Explorer, J. Bohm, AUT 

An Optimization Problem - MS Excel vs CAS, J. Bohm 

Complements de mathematiques (profil Genie Mecanique), M. Beaudin, CAN 

and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
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Volker Loose 


volker@VWLOOSE . DE 


Hello Derive users, 

I entered the following expression 2x + 0.5y + 4z = -2 and performed a multiplication by 3. 
See what happened: 


#1: 2-x + 0.5-y + 4-z = -2 

#2: (2 • x + 0.5-y + 4-z = -2)-3 

3* (4 • x + y + 8-z) 

#3: = -6 

2 

This is correct, but if approximating #2 I get 
#4: 4-x+y+8-z=-4 

Has anyone an explanation for this strange behaviour? 


Joe Frisbee joe.h.frisbee@usa-spaceops.com 

Here is what appears to be happening. Derive 6 does some automatic simplification when using 
’’Approximate”. The second result given is correct but is obtained by multiplying the first result by 
two in order to eliminate the fraction. This can be seen by selecting, from the drop down menus, 
”Options/Mode Settings/Simplification” and then at the bottom of the widow select ’’Display Steps”. 
When ’’Approximate” is used to reduce (2-x + 0.5-y + 4-z = -2 ) • 3 the missing step and its 
explanation will be displayed. I am new to using Derive 6 so I do not know if it is possible to turn off 
this feature. 

#5: Precision := Approximate 

#6: Notation := Decimal 

#7: DisplaySteps := true 

#8: (2 • x + 0.5-y + 4-z = -2)-3 

y 

n-x = y ^ x = 

n 

#9: 4-x+y+8-z=-4 

Andreas Lindner, Austria 

I wanted to demonstrate the influence of parameters in the equation of a plane using slider bars in the 
3D-Plot Window, eg x + 2y+3z = a or ax + by+cz = 4. This does not work. 

DNL. : DERIVE is unable to plot implicitly given expressions in the 3D-Plot Window. Solve the equa- 
tions for z and try once more 


x + 2-y - a 

#3: z — 

3 

a-x + b-y - 4 


#4: 


z 


c 
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Reducing a Matrix - A Chronology of emails 


The Problem: 


Francisco Marcelo Fernandez, Argentina fernande@quimica . unlp . edu . ar 

Dear derivers, 

I have a two-column matrix. In the left column there are some repeated elements that correspond to differ- 
ent elements in the right column. I want to replace the repeated items in the left column by just one element and 
the corresponding different elements in the right by its average. Does anybody have a simple program for doing 
it? 

Example: reduce [1,2; 1,2.1; 1,2. 2; 2, 3; 2, 3. 2] to [ 1 , 2 j- ; 2 , 3 *1 ] 

Regards, 

Marcelo 


John Coleman 


j coleman@GWIA . FRANCISCAN . EDU 


The (hopefully) attached file contains a function which (hopefully) works. It returns the correct answer 
for your example. It is a naive algorithm that makes several passes over the data. A more efficient 
approach is possible of course, if efficiency is an issue. 

Hope that helps 

-John Coleman 

m_reduce_0 1 . dfw 

add_repeats(i , A) 

If A = [] 

#1: 0 

IF(A|1|1 = i, A|l|2, 0) + add_repeats(i , REST(A)) 

count_repeats(i , A) 

If A = [] 

#2: 0 

IF(A|1|1 = i, 1, 0) + count_repeats(i , REST(A)) 

add_repeats(i , A) 

#3: average_repeats(i , A) := i , 

count_repeats(i , A) . 


#4: average_al l_repeats(A) := VECTOR(average_repeats(i , A), i, Ai^l) 

remove_dups(v) := 

If v = [] 

[] 

#5: If MEMBER?(FIRST(v) , REST(v)) 

remove_dups (REST (v) ) 

ADJOIN ( FI RST(v) , remove_dups(REST(v))) 

#6: reduce_matrix(A) := remove_dups(average_al l_repeats(A)) 


1 2.1 


#7: reduce_matri x 1 2.2 


1 2.1 


2 3.1 


2 3 


2 3.2 


to be continued od page 37 






D-N-L#68 


B Kutzler: Optimization Problem with TI-Nspire CAS 


P 5 


A ‘Holistic’ Approach to Solving an Optimization 
Problem with TI-Nspire CAS 


by B Kutzler 


TI-Nspire CAS combines several applications under one roof: A computer algebra system, a graphing 
software, a dynamic geometry system, a spreadsheet, a text editor, and a data & statistics package. 
This makes it an ideal tool for supporting the teaching and learning of topics such as optimization 
problems. Optimization problems can be solved geometrically, graphically, numerically, or algebrai- 
cally. By showing and comparing all four approaches in parallel, students can better understand the 
connections. TI-Nspire CAS is ideal for that. 

We look at the following optimization problem: Given an isosceles triangle with base length ' a ’ and 
height ( h Inscribe a rectangle such that its base lies on the base of the triangle. Among all possible 
rectangles find the one with the largest area. 

We start by drawing an isosceles triangle in the Geometry application. The general Triangle tool is 
not appropriate to begin with, because one cannot enforce that two sides are equally long. Therefore 
we construct it as we would do with pencil, ruler, and compass. We draw a segment, construct its per- 
pendicular bisector, and then choose the segment end points and a point on the perpendicular bisector 
as triangle points. 


3 

1 cm 


ib 



Since the triangle is symmetric with respect to the perpendicular bisector of its base and since the rec- 
tangle’s base should lie on the triangle’s base, the rectangle will have the same symmetry. Therefore, 
again, the Rectangle tool is not appropriate for drawing the rectangle, we need to construct it. We 
choose a point on the triangle base between the left vertex and the base’s midpoint, and then we reflect 
that point at the perpendicular bisector. The point and its reflection make the base of the rectangle. In 
these two points we draw lines perpendicular to the triangle’s base. The intersections of these perpen- 
dicular lines with the two other triangle sides give the other two rectangle vertices. 



Since for the optimization problem we need to measure the area of the rectangle, we need to create the 
rectangle as an object. We use the Polygon tool, which normally would be used for drawing an arbi- 
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trary non-regular polygon. But since the rectangle properties already were enforced, we just use it to 
group the four points into a new geometrical entity. 




We add appropriate labels a , /z, x, and y to the respective parts of the triangle and the rectangle. 



We are ready to measure the lengths of a, h, x, and y, as well as the area of the rectangle. 



Now we have everything we need to find a numeric solution to the optimization problem geometri- 
cally. We grab the left base point of the rectangle, and then move it. We look at the area, and then we 
find a position for the point at which the area is maximal. 




We take a note of all the numbers, or, better, we produce a printout of this page using the File>Print 
command. Next we create a “new” triangle by moving the triangle’s left vertex. 



Again we move the rectangle’s left base vertex to find the rectangle with maximum area for this trian- 
gle. 
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Again we note the resulting numbers. In the classroom we should do more rounds with more triangles, 
and then let the students look at the values of a , /z, x, and y for a pattern. 

Now is a good time to tackle this problem algebraically. We prepare for this by moving the triangle 
into the lower left corner, so that we have more room at the right for a Calculator window. 




We are ready now to solve the optimization problem algebraically. We enter the function of the rec- 
tangle’s area, which is the function to be maximized. We also enter the constraint on the variables 
which stems from the geometric construction. 

area[xy)\=x m y Done 

The constraint on x and y is given by the fact that the upper rectangle points lie on the triangle sides. 
The small white rectangular triangle on the left with sides y and a is similar to the “left half’ of 

the original triangle with corresponding sides h and . Using the theorem of proportional segments 

this gives the following constraint on x and y: y : a ~ x = h : 


y 

_ h 


const!",— 

— 


a-x 

a 


2 

2 



-2 y_ 

2 'h 


x-a 

a 


We solve this equation fory. 

~h u [x-a 

y = — - — 

a 

We turn the solution into an assignment for ‘y’. 


i 

-ki 

i 

ii 

r* 

~h a [x-a) 

a 

a 
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We look at the rectangle’s area, which now should depend only on the variable x (and parameters a 
and h). 


area[x,y) 

i 

'A 

i 

Q 


a 


We differentiate this expression with respect to x, equate the resulting expression to zero, and then 
solve for x. 

. , -h'[2'x-a) 

solve| — -=0,x 


■^-(area(x j/)) 


-#■(2 a x-a) 


a h 

x = — or — =0 
2 a 


We turn the “useful” part of the solution into an assignment for ‘xm’. 

a a 

xm:=— — 

2 2 

We look at y when x takes the value xm. For this we use local substitution by applying the With opera- 
tor. 

y\x=xtn # 

2 

We save this as ym. (We could have combined both steps by entering: ‘y:=ym|x=xm’) 

# h 

ym\=— — 

2 2 

We compute the maximum area. 

area[xm } ym) a u h 

4 

Next we use this formula to check if indeed we found the maximum in the Geometry window. For this 
we make the Geometry window’s measurements available in the Calculator by performing appropriate 
store operations. 


19.0758 cm 2 i 


Store Var 


19.0758 cm 

1 cm 

SlptfElPS 


Unlink 


H := 9.11112 cm 

t 

#=S..j7478 cm 

* 


LinkTo: ► 


A=8. 37478 cm 

A 


If we store this value in a , then a no longer would be a parameter in our TI-Nspire document and the 
calculations in the Calculator would loose (some of) their algebraic nature. Therefore, we should 
change the name to indicate that this is a special value for a. We use ‘aa’. 


19.0758 cm 2 


1 cm 


jj 


:= 9.11112 cm 


#=8,37478 cm 


19.0758 cm 2 

1 cm 

aa=9. 11112 cm 

1 

#=8.37478 cm 



We apply the same procedure to the other three measurements. 
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19,0758 cm 2 ]_ 

aa=9.11112 cm j 
hh=8. 37478 cm / 
XX=4. 56558 
yy=4. 17819 


The current values of aa, hh, xx, and yy now are available in the calculator. 


aa 

9,11112 


aa m hh 

19.0759 

4 



This number is very close to the maximum area which we found in the Geometry window. Need we 
check if our geometric x andy are close to the theoretical results, a ! 2 and h! 2. 


Yes, they are! With our students we may want to change the triangle again, find the maximum rectan- 
gle, and then see how “theory” and “practice”, i.e. algebra and geometry again fit together. 


hh 

4.18739 

2 


yy 

4.17819 


aa 

4.55556 

2 


XX 

4.56558 



More can be done! Let’s take the numbers from the Geometry window and visualize them in a scatter 
plot. The two most interesting numbers for this are ‘xx’, the length of the rectangle base (which is 
controlled via moving the left rectangle base point), and the area. We still have to assign a name to the 
area so that we can use also this value in other applications. We call it ‘ar’. 


We need a Graphing window for a scatter plot - and we need a Lists&Spreadsheet page to collect data 
from a Geometry window and arrange it as a list to be used for the scatter plot. Currently we have one 
page which is split into a Geometry window and a Calculator window. 



ar 

:= 19.0758 cm 2 

1 cm 

aa=9. 11112 cm 



ar=19.Q758 cm 2 

1 cm 

aa=9. 11112 cm 




We want to use two pages for that. The first page should show Geometry and Graphing, the second 
page should show Calculator and Lists&Spreadsheet. This is done easily. We cut the Calculator (with 
( Qt r I ) +( X) after selecting the entire application with the Edit>Select Application command), and 
then add a Graphs&Geometry application which we leave in Graphing View. Then we add a new 
page, split it into two subpages, paste the Calculator into the left half, and then add a List&Spread- 
sheet application in the right half. 


plO 
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£ 

a 

A 

B 

C 

D | 

1 

_ nlvr [-A-(2**-a) 0 


♦ 





solve — u,x 

l a 

| r A =0 

'L 





_ a 


1 





~ 2 

a 

3 





a 

a 

4 





2 

2 

5 





y\x=xm 

h 

| 






2 

1 





Li 

N 

8 




B 

9/99 

^1 



All this is done with only a few mouse clicks! 


Next we prepare the spreadsheet for collecting data from the Geometry window. Column A shall carry 
4 xx’, the length of the rectangle base, and column B shall carry the rectangle’s area. We use the 
Automatic Data Collection command, and then we replace the placeholder ‘var’ by the variable 
for the rectangle base length, which is ‘xx’. 




A 


B 


C 


D 


E 

♦ 

(xx 

► 






We repeat this procedure to collect the rectangle’s area in column B. 


♦ 

A 


=capture(xx. 

=capture(ar,1 

1 

4.56558 

19.0758 




When we move the rectangle vertex again, columns A and B automatically will be filled with the re- 
spective values for ‘xx’ and ‘ar’. But before doing so we will prepare for the values in columns A and 
B to be visualized as a scatter plot. This requires to give names to these two columns so that their con- 
tents will be available in other applications. We choose names ‘base’ and ‘r_area’. (Note that the 
name ‘area’ is already used in the Calculator for the function!) 


♦ 

A 

base 

B r_area 

=capture(xx. 

=capture(ar,1 

1 

4.56558 

19.0758 


In the Graphing window we change the Graph Type to Scatter Plot. This changes the appearance 
of the function editor entry line. 



1 

<3? 


SI 

y*r\ 

A 


From the selection menus we choose base for x and r_area fory. 




Now we have everything together. We make the Geometry window active, grab the rectangle’s left 
base point, and then move it. 
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We change into the Graphing window, and then zoom in on the scatter plot with the Window menu’s 
Zoom - Data command. 



Finally we demonstrate how numbers and algebra fit together. We look at the area function again, . . . 


area[x,y) 

~h m x a {x-a) 


a 


. . . and then produce a specific area function for the current values of a and h (which are available in 
‘aa’ and ‘hh’). This can be achieved using the With-operator. We enter: area(x,y)|a=aa and h=hh 

area[x } })\c!=aa and h=hh 
-.919183"X"(x-9. IIII2) 

The quickest way to plot this function’s graph in the Graphing window is to save this expression as 

fl(x), ... 

/7(x):=-.9 1918281 247769 m* 

Done 

. . . change into the Graphing window, change the Graph Type to Function, in the entry line scroll 
to f 1 , ... 




11 

X 

Cl 

r 

& 




<11111957994) 

X 

A 


. . . and then use the ‘Enter’ -key to plot the graph of fl . 



The essence of this approach is to use the data produced in the Geometry window in all other applica- 
tions and see and study the connections between geometry, algebra, and graphing. 
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Directing Our Suspicions 

Using Matrix Algebra to Organize the Results of an Investigation 

Pat and Carl Leinbach 



Scenario 

You are called by the county police to investigate a shallow grave that has been 
located by two hunters at the edge of the county. You cordon off the scene and 
carefully uncover the remains that are buried in the grave. It is obvious that the 
body was not buried recently. In fact all that remains is essentially bones and 
some shreds of fabric. It is apparent to your experienced eye that the victim 
was a young female, but that is all that you can determine at the scene of the 
grave. You assist in a careful exhumation of the body and direct that it be taken 
to the county morgue. 

Returning to the morgue, you do a more precise analysis. From the length of 
the bones and an examination of the bone ends you determine that the victim is 
a pre-pubescent female. A forensic dental examination reveals that the victim is 
a young girl who was reported missing over 20 years ago. Further analysis re- 
veals that the victim was shot at least three times in the area of the chest. No 
bullets have been found on the scene, but the extent of damage to the bones 
indicates that the weapon was on the order of a .45 caliber pistol. You are un- 
able to make a more exact estimate concerning the weapon. 

At this point your involvement in the investigation is over. Solving the crime is 
the job of the police. A team of detectives is assigned to the case and they be- 
gin the tedious job of interviewing people who knew the victim and also others 
who had contact with her. The town where the victim lived is a small town, so 
the detectives begin there. They also review the local police’s case notes. The 
prime suspect at the time of the disappearance was an African-American 14- 
year old boy who was the son of the girl’s neighbors. His family was the only 
non Caucasian family in the neighborhood. The boy’s sister was the girls “best 
friend”; however, the victim and his sister had just had a small falling out that 
had seemed to have been orchestrated by the victim’s mother. 
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The Police Investigation 

The police began by questioning the now grown boy and his sister. As is the case in so 
many investigations, the questioning revealed other individuals who needed to be visited. 
From their interviews, the police have narrowed the suspect list down to six individuals. Here 
is a very brief summary on the police findings. We will identify the suspects by alphabetical 
letters, A through F. 

A. The African-American boy. Sister, of whom he is very protective, had just re- 
ceived a “rejection” from the victim; Had been seen on the victim’s property the 
night of the disappearance; Determined to have been a fairly normal teen age 
boy, but resentful of the discrimination his family had to endure. 

B. African-American girl. Had been a friend of the victim from the moment her family 
moved next door to the victim’s family amidst a demonstration by the local “citi- 
zen’s committee for a pure neighborhood”; Had recently had a childish falling out 
with the victim; Had written a note to the victim the day of her disappearance ask- 
ing to meet in their “secret place.” 

C. Leader of the citizen’s committee that lead a demonstration greeting the moving 
van when the African-American family moved in; Was very offended by the victim 
being good friend and confidant of suspect B; Quoted as saying “The only thing 
worse than a black is a white person who befriends one”; Known to be a stern fa- 
ther and husband. Known to have kept a gun in his home. 

D. Real Estate agent trying to capitalize on the panic that can result from having a 
non Caucasian family move into the neighborhood; Was trying to buy homes 
cheap and sell houses in his new development on the edge of town; Perceived 
the public friendship of the victim and suspect B as “bad for business.” 

E. Teen aged boy. Son of suspect C; Prejudiced, but confused, especially by the 
friendship of the victim and suspect B; Had been publicly rebuked and slapped 
by the victim when bullying her at school; Apologized to the victim and her friend 
and was, in return, severely disciplined by his father. He was constantly trying to 
please his father. 

F. Mother of victim. Was entranced by suspect D’s attempts to sell her a new 
home, but opposed by her husband; Insisted that her daughter stop seeing sus- 
pect B; Made her daughter ignore her friend in public. 

Each of the above suspects seems to have a reason for being placed on the list of suspects 
to the crime. The job for the detectives is to sort out these reasons and determine on whom 
to concentrate their efforts. 

An old rule of thumb in criminal investigations is to look at three pointers: Motive, Opportu- 
nity, and Profit. Who had the most powerful motive for killing the victim? Who had the best 
opportunity to commit the crime? Who would profit the most from the crime? We will abbre- 
viate these pointers as MOP. In the following investigation, we will apply the MOP pointers 
to the crime at hand. 
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Investigative Objective 

• Organize investigative findings in a sensible way that will lead to the solution of the 
crime. 

Mathematical Objectives 

• Find a way to sensibly incorporate subjective evidence into a meaningful quantita- 
tive analysis. 

• Assign a priority ranking and strength to the various alternatives for solving the 
crime. 

• Using a two way comparison of alternatives to achieve the ranking. 

• Learn about the concept of a dominant vector related to certain matrices. 


Preliminary Mathematical Discussion 

In this section we will introduce the idea of a dominant vector which can be associated with 
certain types of matrices. At the moment, this may seem like an irrelevant mathematical di- 
version, but it really has a lot to do with our attempt to solve the crime; however, the idea of a 
dominant vector (more properly called a dominant eigenvector) is the basis of the Analytic 
Hierarchy Process invented by the mathematician Thomas Saaty in the late 1970’s. This 
process is a way of prioritizing several alternatives into a logical sequence of choices. 

We will illustrate the idea of a dominant vector while working with a 2 x 2 matrix. We will as- 
sume that every entry in the matrix is positive and that the second row is not a multiple of the 
first row. For our example matrix we will use the matrix, 


A = 


7 

1 


8 

3 


We will look at the effect of multiplying the 
matrix consisting of the corners of the follow- 
ing square by the matrix, A (from left hand 
side). 




-1 




-0.5 


-1 

-0.5 

0.5 

1 



--0.5 




■-1 



Since multiplication by the matrix A is a linear process, the lines joining the corners will be 
transformed into lines joining the corners of the images of the corners. Also any point in the 
interior of the square will be transformed into the interior of the image of the square. The 
square can be created by joining the following coordinate points in order: 

Q := [- 72/2, -72/2; -72/2, 72/2; 72/2, 72/2; 72/2, -72/2; -72/2, -72/2] 


We multiply each of the points by A and normalize the resulting vector by dividing it by its 
length. 


D-N-L#68 


Pat & Carl Leinbach: Directing our Suspicions 


pl5 


VECTOR 


f— . V, ,1 

L I A* v | J 


The result is shown as the blue rectangle 
superimposed on the original red square. 
Note that the rectangle is now slanted in a 
particular direction and is a compressed 
form of the original. 


Next multiplying the corner points of the blue 
rectangle by A or doing the same thing by 
multiplying the corner points of the red matrix 
by A 2 , we see a rectangle with a slightly dif- 
ferent slant and even more compressed. 

This rectangle is shown as a green rectan- 
gle. 

This process can repeated as often as you 
like. We show the result after 20 steps multi- 
plying by A 20 and plot the resulting “rectan- 
gle” (black). 





■1 




^ 0.5 N 

* 


\ / 

0 

1 



— 1 



At this stage we have almost reduced the square to a single line. In fact, if we continue tak- 
ing higher powers of A times the coordinates points of the square, we see that the square is 
reduced to a straight line segment lying on the line y = 0.18301x and passing point 
D(0.8453 | 0.1547). 


Exact calculation delivers the coordinates of point D: 






6-2S 2V3-3 
3 ’ 3 


\ 


J 


We choose the vector d (position vector of D ) lying on the line and having the property that 
its components sum to 1 , and call it the dominant vector for the matrix, A. The idea of a 
dominant vector will play an important role in our analysis of the evidence collected by the 
police during their investigation of the suspects in the murder. 


It is important to note that not every matrix has a dominant vector. For example, the matrix 

\/3 _ 1 ~ 

8 = 2 " 2 
1 ^3 
_ 2 2 _ 

has the property that when any vector is multiplied by successive powers of 6 the result is a 
set of 12 different vectors that keeps repeating. 
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73 *x y x -J3 -y 

#9: B-[x, y] = , + 

2 2 2 2 . 

13 ^3 -x y x 73 -y 

#10: B -[x, y] = - , + 

2 2 2 2 . 

k 

#11: VECTOR(B -[2, 2], k, 20) 

The 3 rd expression results in a list of the first 20 products. 

Fortunately for us, we will not be dealing with matrices like B. Instead we will be dealing with 
matrices that have all positive entries. Two German Mathematicians, Oskar Perron (1880- 
1975) and Ferdinand Georg Frobenius (1849-1917), proved that every matrix having all 
positive entries has an associated dominant (Eigen)vector having all positive terms. They 
actually proved much stronger results, but the statement we have given in this paragraph is 
sufficient for what follows. 


Additional comment for students (Josef) 

In the previous paragraph we introduced the expression Eigenvector. This very 
important concept needs some explication. 

With A being a regular quadratic matrix the equation Av = Av can be set up. 
Scalars A and vectors v which fulfil this equation are called Eigenvalues ( A ) with 
their respective Eigenvectors (v). Such Eigenvector only changes its length but 
not its direction if multiplied by matrix A. You can find the Eigenvectors of A using 
DERIVE or the CAS-calculator and then show that one of them is identical with 
the dominant vector from above. 


#47: EIGENVALUES(A) = [2-73 +5,5- 2-73] 

' - 2-73 - 2 ' 

#48: EXACT_EIGENVECTOR(A , 2-73 + 5) = 

-1 

' 2-73 - 2 ' 

#49: EXACT_EIGENVECTOR(A, 5 - 2-73) = 

-1 

[- 2-73 - 2, -1] [ 2-V3 2-V3 

#50: =2 , 1 

I([- 2- 73 - 2, -1]) L 3 3 J 

[- 2- 73 - 2 , - 1 ] 

#51: = [0.84529946, 0.15470053] 

I([- 2-73-2, -1]) 

[2.73 -2,-1] r 2.73 2.73 

#52: = + 2, 1 

I([2.73 - 2, -1]) [3 3 J 

[2-73 - 2, -1] 

#53: = [3.1547005, -2.1547005] 

H[2-73 - 2, -1]) 



The Eigenvectors are presented in 
columns. 


Only this one with equal signs for both 
components can be the dominant Ei- 
genvector. 

Compare the result of #50 with the 
coordinates of point D on page 15. 
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If we apply the comparisons in matrix (1 ) consistently to any ranking of the suspects, then the 
effect of the comparisons should take effect and the rankings should tend towards the line 
determined by the dominant vector. (Try this. Use your CAS to multiply the equal ranking 
vector (1,1, 1,1, 1,1) by the 30-th power of the Motive matrix and normalize the result by divid- 
ing by the sum of the terms so that they add up to one). If you have down loaded the ac- 
companying DERIVED file called , AHP_Analysis . dfw, the program, Dom_vec, will find the 
dominant vector for any positive matrix up to a given tolerance between the actual and com- 
puted value. Here is the result for the Motive matrix. Compare it with your result. 

#6: Dom_vec (Motive, 0.0001) 

#7: [0.129723, 0.145604, 0.0868161, 0.12758, 0.418666, 0.0916085] 

Suspect E, corresponding to the fifth value in the list is clearly the highest ranked according 
to this criterion. Suspect B is a distant second followed by suspects A, D, F, and C. 

Additional comment for students (Josef) 

We come back to the Eigenvector and show that our approximative procedure 
can be confirmed in a very mathematical way by applying this concept. 

#8: EIGENVALUES(Motive) 

#9: [9.52509, 0.144898, -1.85023 + 4.74921 -L, -1.85023 - 4.74921- i, 0.0152352 

+ 2 . 87102 • i , 0.0152352 - 2. 87102- i.] 

#10: dv := APPROX_EIGENVECTOR(Motive , 9.52509) 

#11: dv := [-0.261811, -0.293876, -0.175191, -0.257476, -0.844939, -0.184882] 

dv 

# 12 : 

I(dv) 

#13: [0.129726, 0.145614, 0.0868066, 0.127578, 0.418664, 0.0916085] 

DERIVE returns - approximatively - two real Eigenvalues (approximatively be- 
cause the characteristical equation is of degree 6). The attempt with the first real 
Eigenvalue leads to the appropriate Eigenvector. In case of larger matrices calcu- 
lating the Eigenvalues can lead to problems. It seems to be better to use an itera- 
tive process to approach the “limit”. Beside the offered program you can perform 
the first, say 20, multiplications and either observe the whole series of vectors or 
just the last row of the resulting table: 

#13: [0.129726, 0.145614, 0.0868066, 0.127578, 0.418664, 0.0916085] 

f ( Motive -u 

ITERATES , u, [1, 1, 1, 1, 1, 1], 20 

#14: v l I(Motive-u) 

21 

It is no problem to find the Eigenvectors - even for larger matrices - on the CAS- 
calculators without calculating the Eigenvalues. 
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The four detectives compare the suspects wirh respect to the remaining criteria which 
are Opportunity and Profit and calculate the respective dominant vectors. 



#18: Dom_vec (Opportunity , 0,00001) 

#19: [0,133494, 0.287197, 0,0399400, 0,0183003, 0,0904892, 0.430578] 

#20: Dom_vec(Prof it , 0.00001) 

#21: [0,0223754, 0,0667261, 0,151181, 0.486502, 0,212729, 0,0604838] 

Confirm the results by calculating the Eigenvectors. 

Upon examining the dominant vectors for the three rankings, you see that each criterion 
gives a different ranking of the three suspects. The team of detectives has one more “collec- 
tive decision” to make: How does the evidence fit the crime they are investigating? In other 
words, how much weight should they give to each of the three criteria relative to the circum- 
stances of this crime? Once again, they can use the Analytical Hierarchical Ranking 
Scheme to assign the weights. 
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' 1 5 8' 

1 

14 

#22: Priorities := 5 

1 1 

.84 

Try to explain the elements of this matrix. 

With Dominant Vector: 

#23: Dom_vec(Pri oriti es , 0.00001) = [0.733384, 0.199069, 0.0675455] 

It is clear that Motive and Opportunity are the dominant criteria in this case and the Profit has 
relatively little role to play in the opinion of the detectives. This should be expected in the 
murder of a child coming from a middle class family. Indeed, Motive is by far the prime crite- 
rion. Each individuals score for the three criteria is multiplied by the weighting given to the 
criteria to find a final score for each suspect. 

' 0.129727 0.133494 0.0223754 1 [ 0.123225 ' 

0.145614 0.287197 0.066726 0.16847 

[ 0.733384 1 I 

0.0868068 0.03994 0.151181 0.0818251 

#24: • 0.199069 = 

0.127578 0.0183003 0.486502 0.130067 

L 0.0675455 J 

0.418664 0.0904892 0.212729 0.339423 

. 0.0916082 0.430578 0.0604837 J L 0.156984 . 

The rows of the left matrix are formed by the values for Motive, Opportunity and Profit for the 
persons A though F. The product delivers the final evaluation, which must be sorted. As the 
5 th row is showing the maximum, the teen aged boy seems to be the main suspect. At the 
other hand the minimum value in the 3 rd puts the committee leader - which is suspect C - on 
the last position of suspects. 

This can be done using the DERIVE E® program MOP_Suspects that is a part of the package 
supplied with this activity. The result of the program is given below. 

' 0.339425 E ' 
0.168469 B 
0.156985 F 

#25: MOP_Suspects (Priori ties, Motive, Opportunity, Profit, 0.00001) = 

0.130068 D 
0.123225 A 
. 0.0818258 C . 

Upon further investigation and questioning, it turned out that suspect E had committed the 
murder and confessed to the murder. He was arrested, tried, and found guilty. 
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Exercises Related to the Activity 

1. Read a mystery novel and prior to reaching the end of the novel, identify the ma- 
jor suspects. Use the Analytical Hierarchical Ranking Method to identify the pos- 
sible perpetrator of the crime. Compare your result with the author’s solution of 
the crime. 

2. Together with a small group of friends, watch a detective show on television that 
has a good development of the evidence. Have everyone leave before the end of 
the show, but videotape the ending portion that you missed. You and your 
friends create the matrices for the Analytical Hierarchical Ranking Method and 
determine the most likely guilty suspect. Watch the end of the show and com- 
pare your result. 

3. The Analytical Hierarchical Ranking Method has applications in many areas. Ap- 
ply the method to making an important decision in your life, such as choosing a 
college to attend, choosing a summer job, choosing a destination for a family trip, 
etc. Obviously you will use criteria other than the ones used for crime solving. 

4. Choose important criteria for success in a particular sport. Rank the teams in 
your league relative to these criteria using the Analytical Hierarchical Ranking 
Method. How do the results compare to the final rankings of the team after the 
end of a season? 

Final Comments of the Coroner 

The role of the coroner ended in this case at the point the coroner assigned an manner 
(murder) and means ( gun shot wound to the chest) to the death of the victim. The remain- 
der of the work was for the police. They did the interviews and gathering of evidence. 

Crime solving, like the practice of medicine, is not a science. Scientific methods are applied 
to the evidence and conclusions are reached. However, it is seldom that the results are 
based on scientific laws or indisputable facts. There is a lot that requires the insight and 
training of the investigator. The methods used in this activity are a way for the investigators 
to organize their insights and “hunches”. They are nothing more. They can be used to direct 
an investigation, not produce irrefutable evidence. In the case of the above investigation, 
this organization produced a focus that lead to a confession by the guilty party. 


More information about AHP can be found on the following websites: 

http://en.wikipedia.org/wiki/Thomas_Saaty 
http://en.wikipedia.org/wiki/Analytic_Hierarchy_Process 
http://www.boku.ac.at/mi/ahp/ahp.pdf (deutsch) 

http://www.boku.ac.at/mi/ahp/ahptutorial.pdf (english) 
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Groebner6.dfw - a DfW Utility File for computing Grobner Bases in 

DERIVE 6 

Alessandro Perotti, Dept, of Mathematics, University of Trento, Italy 
perotti@science.unitn.it 
www.science.unitn.it/~perotti/groebner.htm 


This DfW utility file contains a function that computes Grobner bases of a set of polynomials 
with respect to a monomial ordering and other related functions. DERIVE 6 contains the new 
internal function GROEBNER_BASIS to construct a Grobner basis for a collection of polyno- 
mials based on lexicographic ordering of the variables. We introduce new functionalities giv- 
ing the possibility of choosing the monomial ordering. Besides allowing the explicit com- 
putation of Grobner bases and normal forms, the new functions can be used, for example, to 
eliminate some variables between the equations and to perform the Gaussian elimination of 
a system depending on parameters. 

We wish to thank V. Anisiu for his useful suggestions. 

• 0. Introduction 

• 

A Grobner basis is a set of polynomials in one or more variables that has special properties, 
useful to solve many problems involving polynomial equations. This concept generalizes 
Gaussian elimination for linear equations and the Euclidean algorithm for computing the 
GCD of two polynomials in one variable. Using the words of B.Sturmfels [2], "Given that 
polynomial models are ubiquitous across the sciences and engineering, Grobner bases have 
been used by researchers in optimization, coding, robotics, control theory, statistics, molecu- 
lar biology, and many other fields." 

To give some examples of Buchberger's Algorithm for the construction of a Grobner basis, 
we start with a system of two linear equations in three variables 

3x + 2y + z = 2, 2x - y - 4z = 1 . 

We rewrite it as f| = 3x + 2y + z - 2 = 0, f 2 = x - y - 4z - 1 = 0. To eliminate the variable x 
from the equations, we form the new linear polynomial 

f 3 = 1/3 U - 1/2 f 2 = 7/6 y + 7/3 z - 1/6 

and get the equivalent set of equations ^=0, f 2 =0, f 3 =0. What we have done is: order the 
variables as x>y>z, then use the leading terms 3x of fi and 2x of f 2 to eliminate the first vari- 
able x and get the equation f 3 =0. 

The set {f^ f 2 , f 3 } that we have obtained is a Grobner basis of the initial set {f-i, f 2 }, since the li- 
near combinations of the leading terms of f-i, f 2 , f 3 give all the possible combinations of the 
leading terms of the original polynomials fi, f 2 . From the Grobner basis {fi, f 2 , f 3 }, a reduced 
Grobner basis can be obtained. In the linear case it corresponds to the row-reduced echelon 
form of the system of equations. In our example a reduced basis is given by 


{gi = x - z - 4/7, g 2 = y + 2z - 1/7} 


p24 


Alessandro Perotti: Grobner Bases in DERIVE 6 


D-N-L#68 


In the non-linear case, in order to make a similar construction, we have to fix an ordering for 
the monomials. Given an ordering for the variables, there are many orders which are compa- 
tible with the product of monomials. For example, the lexicographic ordering (used by the in- 
ternal function GROEBNER_BASIS) for monomials in two variables x > y is 

1 < y < y 2 < y 3 < ... < x < xy < xy 2 < xy 3 < ... < x 2 < x 2 y < ... < x 3 < ... 

while the graded lexicographic ordering gives 

1 < x < y < x 2 < xy < y 2 < x 3 < x 2 y < y 3 < ... 

Once the monomial ordering is fixed, every polynomial has a leading term (a monomial mul- 
tiplied by a coefficient). For example, let F be the set 

{fi = x 2 + y 2 + z 2 + w 2 , f 2 = x + z 3 - w 3 , f 3 = x 2 + 2 y 2 - y z - w 2 }. 

The leading terms of f|,f 2 ,f 3 w.r.t. the graded lexicographic ordering are, respectively, x 2 , z 3 , 
x 2 . Now we take elements of F in pairs and cancel out their leading terms by means of poly- 
nomial combinations of them. We form the so-called S-polynomials 

U = z 3 h - x 2 f 2 = - x 3 + w 3 x 2 + z 3 (y 2 + z 2 + w 2 ) 

f 5 = fi - f 3 = - y 2 + y z + z 2 + 2 w 2 

f 6 = x 2 f 2 - z 3 f 3 = x 3 - w 3 x 2 - z 3 (2 y 2 - y z - w 2 ) 

We now apply the division algorithm and find that 

U = (w 3 - x) fi + (y 2 + w 2 + z 2 ) f 2 

and then we don't need to add f 4 to the list {fi, f 2 , f 3 }. Instead, the leading term of f 5 , which is 
-y 2 , can not be cancelled out by means of the leading terms of fi, f 2 , f 3 .Then we add f 5 to the 
set F and get F' = {fi, f 2 , f 3 , fs}. When we apply the division algorithm to f 6 and the new set F' 
we get 

f 6 = (x - w 3 ) fi + (- y z - 2 z 2 - 3 w 2 ) f 2 + (x + 2 z 3 - w 3 ) f 5 
and so we don't add f 6 to the set F. 

If we now form the S-polynomials of pairs containing f 5 and divide them through the elements 
of F, we get only polynomial combinations of fi, f 2 , f 3 , f 5 . 

Buchberger's Criterion says that F is then a Grobner basis of the set F the polynomial com- 
binations of the elements of F are the same obtained by F and the leading terms 

x 2 , z 3 , -y 2 

of the elements of F divide the leading terms of every polynomial combinations of the ele- 
ments of F. 

From F the unique reduced Grobner basis G of the set F can be obtained applying the divi- 
sion algorithm. It is a Grobner basis characterized by two properties: every leading term has 
coefficient 1 and every non-leading term of the elements of G is not divisible by the leading 
monomials of F (which are the same as those of G). 

In our example, the reduced Grobner basis w.r.t. the graded lexicographic ordering is 
G = {gi = y 2 - y z - z 2 - 2 w 2 , g 2 = x 2 + y z + 2 z 2 + 3 w 2 , g 3 = x + z 3 - w 3 } 
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with leading terms y 2 , x 2 , z 3 . The Grobner basis reveals geometric properties of the set V of 
common zeroes of The absence of a power of w in the leading monomials of G means 
that there are infinite common zeroes (over the complex numbers). More precisely, it says 
that V is one-dimensional. This is similar to the linear case, where the number of not-leading 
variables in the reduced echelon equations gives the dimension of the linear variety of ze- 
roes of a compatible system. Besides, the set V is empty exactly when the reduced Grobner 
basis is {1}. 

Note that the internal function GROEBNER_BASIS returns a Grobner basis w.r.t. the lexico- 
graphic ordering which is not necessarily reduced, since the coefficients of the leading terms 
can be different from 1 . 


• 1. Functions 

Groebner ( f , v , ord) 

Eliminate ( f , vl , v2 ) 
TotalDegree ( f , v) 

OrdMatrix (ord, n) 

E10rd(k, n) 

LeadingTerm ( f , v, ord) 

NF (f , ff, v, ord) 

Buchberger (f , v, ord) 

System (e, v, k) 

Eigensystem(a, ord, v, w) 
RREF (a , k) 

RREF (a , b, k) 

LagrangeMultipliers (f ,g,v) 


reduced Grobner basis for a set of polynomials 
elimination of a set of variables in equations f=0 
total degree of a polynomial 
matrices inducing relevant monomial orderings 
matrices inducing elimination orderings 
leading term of a polynomial 

normal form of a polynomial with respect to a 
set of polynomials 

(non reduced) Grobner basis for a set of 
polynomials 

gaussian reduction of a system depending on 
parameters 

reduction of the eigensystem Ax=xx 
row-reduction of a matrix depending on 
parameters 

row-reduction of augmented matrix depending 
on parameters 

computation of min-max points using Lagrange 
multipliers 


• 2. Description of the functions 

• 2.1 Computation of Grobner Bases 

Groebner ( [fl , . . . , fm] , [xl , . . . ,xn] , ord) 
computes the reduced Grobner basis for the set of polynomials generated by f 1 , . . . , fm 
with respect to the ordered variables xl, . . . ,xn and to the monomial ordering described 
by ord. 
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The argument ord that defines the monomial ordering can be: 

- an integer between 1 and 4 (or the variables lex, grlex, grevlex, invlex), correspond- 
ing to lexicographic, graded lexicographic, graded reverse lexico- 
graphic, inverse lexicographic 

- a non-singular matrix with non-negative integral entries that induces the ordering (see [1] and 2.2 for 
more information). 

• The default ordering is lexicographic with respect to the variables xl>x2> . . . >xn . 

• If also the second argument is omitted, the ordering is lexicographic with respect to the 
variables appearing in f 1 , . . . , fm ordered by internal Derive ordering. 


• Examples 


[22222 2 2 3 3] 

#1: Groebner(Lx +y + z +w,x + 2-y - y-z -w,x+z -wj, 


[x, y, z, w], grevlex) 


#2: 

#3: 

#4: 


[2 222 22 3 3 1 
Ly - y-z - z - 2-w , x + y-z + 2-z + 3-w , x + z - w J 

[22222 2 2 33 1 

Groebner(Lx + y + z + w,x + 2-y - y-z - w,x + z - w , w - lj) 


12 9 8 6 5 4 3 2 

w - 1, z -4-z +5-z +12-Z -10-z +5-z -16-z +18-Z +16, 


#5: 


#6: 


11 7 4 3 

z 8 5 • z 5 5 • z 5-z 2 5-z 3 

y - + z - - 2-z + - + 2-z - , x + z - 1 

4 4 2 4 2 

[22222 2 2 33 1 

GR0EBNER_BASIS(|_x + y + z + w,x + 2-y - y-z - w,x + z - w , w - 1 J , 

[x, y, z, w]) 


[w - 1, 


12 98 6 54 3 2 

z -4-z +5-z +12-Z -10-z +5-z -16-z +18-Z +16, 


11 8 7 5 4 3 2 3 1 

4-y - z + 4-z - 5-z - 8-z + 10-z - 5-z + 8-z - 10-z, x + z - lj 


#7: 


#8: 


f\ 2 3 2 4 1 1 3 T| 

Groebner Lx -y - x-y , x -y - xj, [x, y], 

l L 3 0 JJ 

f 3 2 1 

Lx - x, x-y - x J 


Eliminate ( [fl fm] , [xl,...,xh], [yl,...,yk]) 

constructs equations in which some variables have been eliminated. It eliminates the first set 
of variables xl, . . . ,xh between the polynomial equations fl=0, . . . ,fm=0 in the vari- 
ables xl , . . . , xh , yl , . . . , yk . 
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• If the third argument is omitted, the equations are considered with respect to all the variables 
appearing in f 1 , . . . , fm. 

• As it is shown in the examples of section 3, the function can be used in some cases also when the 
equations are not polynomial. 


• Examples 



1 

r 2 2 

2 2 2 2 2 3 3 ] 

#9: 

Eliminated 

Lx + y 

+ z + w , x + 2*y - y*z -w,x+z - w , w - 1 J , [x, w]) 


r 

2 

2 6 3 2 1 

#10: 

L: 

/ - y*z 

- z - 2, y*z + z - 2*z + 2*z + 4j 


1 

r 2 3 

2 4 1 

#11: 

El imi nate(| 

Lx -y - 

X 

X 

X 

1 

X 

1 — 1 
X 

#12: 



[] 


• 2.2 Other functions related to Grobner Bases 

TotalDegree (f , [xl,...,xn]) 

computes the total degree of a multivariate polynomial f in the variables xl , . . . ,xn. 

5 7 5 

#13: Total Degree(2-x -y - 3-x-y-z ) 

#14 : 12 

OrdMatrix (ord, n) 

gives an integer square matrix of order n that induces the rational monomial ordering ord, 
which can be lex or 1 , grlex or 2, grevlex or 3, invlex or 4. See [1] for more 
details. 

#15: MAP_LIST(OrdMatrix(ord, 3), ord, [1, 4]) 



■ 1 

0 

i 

O 


' l 

l 

1 

O 


' 1 

1 

1 ‘ 


1 

O 

0 

I— 1 

1 



0 

1 

0 

j 

l 

0 

l 

j 

1 

1 

0 

j 

0 

1 

0 



1 

o 

0 

l . 


. l 

0 

— i 

o 


1 

0 

1 

O 


. l 

0 

1 

O 



E10rd(k, n) 

gives an integer square matrix of order n that induces an elimination ordering with respect to 
which the first k variables always precede the others. 

#17: E10rd(2 , 4) 


#18: 


'10 0 1 ' 
10 0 0 
0 110 


0 


1 


0 


0 
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LeadingTerm(f , [xl,...,xn], ord) 

computes the leading term of the polynomial f with respect to the ordering ord in the vari- 
ables xl>x2>. . . >xn. 


• The default ordering is lexicographic with respect to all the variables appearing in f . 


2 3 2 2 2 5 

#19: MAP_LIST(LeadingTerm(2-x -y + 3-x -y -z - x-y , [x, y, z], ord), ord, 

[lex, grlex, grevlex, invlex]) 


#20: 


23 222 5 22 21 

2-x -y , 3-x -y -z , - x-y , 3-x -y -z J 


NF(f , [fl, . . . ,fk] , [xl, . . .xn] , ord) 
returns the remainder of the polynomial f on division by f 1 , . . . , f k with respect to the or- 
dering ord. 

If [fl fk] is a Grobner basis, the remainder does not depend on the ordering of the polynomi- 

als and is called Normal Form of the polynomial f with respect to [f 1 , . . . , fk] . It is zero if and only 
if f is a polynomial combination of f 1 , . . . , f k . 


• The default ordering is lexicographic with respect to all the variables appearing in f . 


#21: 

NF(3-x 2 -y, 

[ 3 3 2 2 2 

Lx + z - w , x + y*z + 2*z + 3*w , 

2 2 2] 

, y - y-z - z - 2 • w J, 


l 1 

X 

N 

w] , grlex) 


#22: 


2 2 2 

3*x - 9«y«(z + w ) - 6«w *z - 3- 

3 
• w 

#23: 

NF(3-x 2 -y, 

[ 3 3 2 2 2 

Lx + z - w , x + y*z + 2*z + 3«w , 

2 2 2] 

, y - y-z - z - 2 • w J) 

#24: 


6 3 3 6 

y • (3 *z - 6*w »z + 3 *w ) 



Buchberger ( [fl , . . . , fm] , [xl , . . . ,xn] , ord) 
implements the Buchberger algorithm for the computation of a (non-reduced) Grobner basis 
for the set of polynomials generated by f 1 , . . . , fm with respect to the ordered variables 
xl, . . . ,xn and to the monomial ordering described by ord. The basis obtained always 
containes the generators f 1 , . . . , fm. 


• The function Groebner (f ,v,ord) uses an improved version of Buchberger algorithm for the 
computation of the reduced basis. 

• The default ordering is lexicographic with respect to all the variables appearing in f . 

[22222 2 2 3 3] 

#25: Buchberger([x +y +z +w,x + 2-y - y-z -w,x+z -wj, [x, y, z, w] 

[x, y, z, w], grevlex) 


[22222 2 2 332 2 2] 

#26: Lx +y +z +w,x +2-y -y-z-w,x+z - w , - y +y-z+z +2-wJ 

[22222 2 2 3 3] 

#27: Buchberger([x +y +z +w,x + 2-y - y-z -w,x+z -wj) 
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222222 2 3 3 2 22 


+ y + z 

+ w 

, x + 2-y - y-z - w 

, X + z - w 

> - y 

+ y-z + z 

+ 2-w , 

6 

3 

3 2 6 2 

2 4 


11 3 

8 7 

y-z + z - 

2-w • 

z + 2-z + w + 3-w , 

- w -y • (w + 

3) + 

z - 4-w - 

•z + 5-z + 

5 6 

2 

3 4 3 

5 2 4 


6 

2 

z -(5-w + 

3 • w ) 

- 10-w -z + 5-z - 2 

•w -z • (w + 

3) + z 

-(3-w + 7 < 

•W ), 


17 3 14 13 11 6 2 3 10 9 5 8 4 

z - 6-w -z + 5-z + z -(14-w + 6-w ) - 20-w -z + 5-z - 8-w -z -(2-w + 

7 6 2 365 12 8 4 544 

3) + z -(25-w + 13-w ) - 10-w -z + z -(9-w + 30-w + 9-w ) - 2-w -z -(5-w + 

7 2 8 4 14 3 11 10 2 8 4 

+ 13) - 2-w -z • (w + 6-w + 9), - z + 4-w -z - 5-z - 6-w -z • (w + 1) + 

37 659 5244 428 4 

10-w -z - 5-z + z -(4-w + 12-w ) - w -z -(5-w + 13) - w -z • (w + 6-w + 9), 

3 13 12 6 10 3 9 8 7 9 5 

2-w -z - 2-z - 8-w -z + 18-w -z - 10-z + z -(12-w + 12-w ) - 

264 354 12 8 39 5 

4-w -z -(8-w + 3) + 30-w -z - 2-z -(4-w + 12-w + 5) + z -(18-w + 50-w ) - 

2 2 4 15 11 7 12 8 4 

2-w -z -(5-w + 13) + z-(2-w + 12-w + 18-w ) - 2-w - 12-w - 18-w , 

12 39 866 2 35 4 534 

z - 4-w -z + 5-z + z -(6-w + 6-w ) - 10-w -z + 5-z - 4-w -z • (w + 3) + 

2 6 2 12 8 4] 

z -(5-w + 13-w ) + w + 6-w + 9-w J 


• 2.3 Elementary applications to linear systems and matrices 

System( [el , . . . ,em] , [xl , . . . ,xn] , [kl , . . . ,ks] ) 
performs the gaussian reduction of the system of linear equations el=0, . . . ,em=0 in the 
variables xl , . . . ,xn depending polynomially on parameters kl , . . . ,ks . 

If the vector of parameters is omitted, then in the resulting equations may appear rational functions of 
the parameters and they may be not equivalent to the given equations for some values of the parame- 
ters. 


#29: System(|_x + y + z - 1, k-x + 4-y + 9-z + 11, k-x + 2-y + 3-z + lj , [x, y, z], 

[k]) 

#30: [z-(k - 3) + 3-k - 9, y-(k - 2) - 4-k + 8, y-(z + 3) - 4-z - 12, x + y + z - 1] 

From the reduced system we get that if k=3 the system has infinite solutions x = -3-t, y = 4, z 
= t (t any real), if k=2 the system has infinite solutions x = 4-t, y = t, z = -3 (t any real), else it 
has a unique solution x = 0, y = 4, z = -3. 

#31: System([x + y + z - 1, k-x+ 4-y + 9-z + 11, k-x + 2-y + 3-z + l], [x, y, z]) 


#32: 


[z + 3, y - 4, x] 
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Eigensystem (a, ord, [xl, . . . ,xn] , w) 
performs the gaussian reduction of the eigenvalues system Ax=Xx associated to a square 
matrix a. The default ordering is lex, with default linear variables xl, . . . ,xn and eigen- 
value variable w. 


f 2 l o 


#33: Eigensystem 


13 4 


L 0 0 2 


#34: Lw*x3 - 2*x3, w *x2 - 5*w*x2 + 5-x2, x2*x3, - w*x2 + xl + 3*x2 + 4-x3 



s 

' 2 

1 

0 ' 

\ 

#35: Eigensystem 


1 

3 

4 

, grlex, [x, y, z], t 


V 

. 0 

0 

2 . 

/ 


#36: Lz.(t - 2), -x + y-(t - 3) 


2 

4-z, y-z, x • (t - 2) - y, x-z + 4-z , 


2 2 
x + x-y - y 


2 

16 -z 


] 


RREF (a , [kl , . . . ,ks] ) 

tries to compute a Row Reduced Echelon Form of a matrix a depending on parameters 
kl , . . . ,ks. 

RREF (a , b , [kl f . . . ,ks] ) 

computes the Row Reduced Echelon Form of the augmented matrix [a,b] depending on 
parameters kl , . . . ,ks. The second argument b can be a vector or a matrix. Any null row is 
deleted. 

If the parameters are omitted, then the internal function row_reduce is called for a more efficient 
computation. In this case, the reduced matrix may contain rational functions of the parameters and the 
reduction may be not valid for some values of the parameters. 



r 

'2 3 4-h 

N 

#37: RREF 


3 2-h 1 

, [h] 


K 

.02 1 . 

> 


#38: 


5 

1 0 

28 

1 

0 1 

2 

13 

0 0 h 

28 . 


The RREF shows that the matrix has rank 2 for h=1 3/28, rank 3 otherwise. 


2 3 4 • h 

3 2-h 1 

0 2 1 


, [2, 1, 1], [h] 


#39: RREF 
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#40: 


5 9 2 • h 

1 0 

28 28 7 

1 1 

0 1 

2 2 


13 h 1 

0 0 h 

28 7 28 


2 3 4-h 1 


#41: RREF 


3 2 • h 1 

0 2 1 . 


#42: 


'10 0 ' 
0 10 
. 001 . 


• 2.4 Another application: the Lagrange multipliers 


LagrangeMultipliers (f , [gl, . . . ,gk] , [xl, . . . ,xn] ) 
applies the method of Lagrange multipliers for the computation of critical points of a polyno- 
mial function f (xl, . . .xn) subject to k polynomial constraints gl=0, . . . ,gk=0. It returns 
a system of equations reduced with respect to lexicographic ordering, whose solutions are 
the coordinates of the critical points. 


2 

#43: LagrangeMul ti pi iers(z 


x*y*z + x 




Y-z 


-2L 


#44: 


7 4 

8 6 z 5 z 2 

z - 4-z - 4, y - + 2-z , x - + 2-z 

2 2 


2 

#45: LagrangeMul ti pi iers(z - x-y-z + x, [y-z - 2]) 


#46: 


[ 1 ] 


If the constant polynomial 1 is returned, the function has no critical points subject to the con- 
straints. 


• 3 Other examples 

In the univariate case, the Grobner basis contains only the GCD of the polynomials. 

[9 87 6 5 4 3 2 

#47: Groebner(Lx - 3-x + x - 3-x - 3-x + 6-x + 17-x - 22-x - 11-x + 15, 

7 6 5 4 3 2 1 

3-x - 9-x + 5-x - 15-x - 4-x + 3-x + 48-x - 63 J) 


#48: 


[x - 3] 
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3 2 64223 2 543 2 

z - 2*z +2)+y -y*z +y • (z - 2*z + 1) - z -z +z + 3*z - 2, 

2432 65432 

x«Cy *(3*z -z -6*z +4)+z - 3*z - 3*z +3*z +6*z -4)+ 

86432 24 2 64 

y - y + y *(2*z - 5-z + 3) - y -(7-z - 14-z + 7) - 2-z + 9-z - 

3 2 1 

2-z - 9-z + 4 J 

An assignment to a variable is a convenient way for using a lengthy result in other computa- 
tions: 

#61: PROG(gb := Groebner(polys , [x, y, z] , lex), MAP_LIST(LeadingTerm(f) , f, gb)) 

[ 12 11 2 4 2 4 2 2 31 

#62: Ly , x-z , 24-x-y -z, x-y , 12*x *z , x *y , x J 

This means that the lex Grobner basis contains 7 polynomials... 

#63: MAP_LIST(DIM(TERMS(EXPAND(f))) , f, gb) 

#64: [25, 49, 53, 9, 49, 6, 4] 

...and that these polynomials contain up to 53 terms... 

#65: MAP_LIST(TotalDegree(f) , f, gb) 

#66: [12, 13, 12, 5, 12, 4, 3] 

...and their maximum total degree is 13. 

10 3 3 

#67: NF(x -y -z , gb) 

11 37 6 3 39 63 

#68: y -z + y -(2-z -2-z)+y-(z -2-z +z) 
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I'd like to add two examples demonstrating the power of Groebner Bases (Josef) 

In DNL#25 Karl Heinz Keunecke presented a system of linear equations 

[2-x + 6-y - 3-z +6=0, 4-x + 3-y + 3-z -6=0, 4-x - 3-y + 9-z - k = 0] 

which - according to the result of applying ROW_REDUCE - seems to have no solution: 


2 6-3-6 


R0W_REDUCE 4336 


4 -3 9 k 


0 1-1 0 


0 0 0 1 


The Groebner Basis shows that there is a solution: 

GR0EBNER_BASIS([2 • x + 6-y - 3-z + 6, 4-x + 3-y + 3-z - 6, 4-x - 3-y + 9-z - k] , [x, y, z, k]) 

[k - 18, y-z+2, 2-x+3-z-6] 

For k = 18 we have x = ^ , y = z - 2. 

2 ' 

Let's try to solve a polynomial system by applying Groebner Bases. The following system is to be 
solved: 

r 2 2 2 2 2 2 1 

equs := Lx + y =5, x + z + y = 10, y + x -z = 20 J 


We use the implemented GROEBNER BASIS function and we get 

[ 2 2 2 2 2 2 1 
grb := GROEBNER_BASIS(Lx + y - 5, x + z + y - 10, y + x -z - 20 J, [x, y, z]) 


[6 5 4 

grb := Lz - 20 -z + 93 -z + 7C 


4 3 2 2 

+ 605, 15 -y - z + 20 - z - 94-z - 35-z + 230, 15-x 


+ z - 20 - z + 94 


2 

-z + 50 - z - 330. 


The first equation can be solved for z: 

zval5 := SOLUTIONSCgrb , z) 

1 

zvals := [7.557203220, 12.10993346, 1.634943157 + 0 . 1303137913 - L , 1.634943157 - 0 . 1303137913 - l , 
-1.513511499 + 0. 09323222015 -l, -1.513511499 - 0 . 09323222015 - L ] 

We substitute the 6 solutions for z in the 2 nd equation which then is solved fory: 

yvals := VECTOR (SOLUTIONS (v , y), v, VECTORfgrb , z, zvals)) 1 

2 


yvals := [2.175471042, -2.212919732, 0.5373527040 + 1. 751123793- i, 0.5373527040 - 1 . 751123793 - i , 
0.4313716390 - 2 . 503763362 - L , 0.4313716390 + 2 . 503763362 - i ] 


We substitute the 6 solutions for z in the 2 nd equation which then is solved for x: 
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Surface #3 : x 2 + y 2 = z 2 - z 3 

2 2 2 3 

VECTOR(Contou rPts_XY (x + y - z + z , 1, -1.5, 1.5, -1.5, 1.5, 0.05, 0.05), 1, -1, 1, 0.1) 



2 2 2 3 

Imp! i ci tDots(x + y = z - z , [-1.5, -1.5, -1], [1.5, 1.5, 1], 0.05) 



Using polycontour . mth (DNL#63) and IMPLICIT_Peter.mth (DNL#64) 



Plot with DPGraph. The right figure shows the plot of ax 2 + by 2 = c(z 2 - z 3 ) varying the parameters a , 
b and c using scroll bars. 
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FAZ - Flache #3 



Autograph 3.20 


MuPad 


User Form (continued from page 4) 

John Coleman 

Greetings, 

The function I came up with yesterday performs (as I suspected) poorly on large examples. For exam- 
ple, on my machine it took 124 seconds to evaluate reduce_matrix(A) when A was given by A:= vec- 
tor([random(10),random(10)],i, 1,1000). The attached file contains a programmed version which takes 
less than a second with the same example. Please let me know if you encounter any bugs. 


m reduce 02.dfw 


# 1 : 


pre_process(v, u, f, i, j, n, e, x, w) 

Prog 
u := [] 
f := [] 
i 1 
n := DIM(v) 

Loop 

If i > n 
RETURN u 
e := vfi 
x efl 

j :z POSITION (x, f) 

If j = false 
Prog 

f APPEND(f, [x]) 

u APPENDCu, [ [x , ef2, 1]]) 

Prog 
w := 14 j 

w |2 w |2 + ej ,2 

Wf3 Wf3 + 1 

14 j w 
i := i +1 


av_al 1 (A) VECTOR(av(A ), i, 1, DIM(A)) 
#3: i 

#4: reduce_matrix(A) av_al 1 (pre_process(A)) 

#5: A VECT0R([ RANDOM (10) , RANDOM(IO) ] , i, 1, 

1000) 


# 6 : reduce_matrix(A) 


3 4.59047619 

7 4.551401869 

2 4.972222222 

4 4.885714285 

9 4.842105263 

1 4.740384615 

8 4.420454545 


# 2 : av(u) 



6 4.034090909 
0 4.051020408 
5 4.705882352 


0.141 sec 


-John Coleman 
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Johann Wiesenbouer j . wiesenbauer@TUWIEN. AC .AT 

Hello John, 

Yes, it is true, your second version is considerably better than the first one. Even though its 
performance breaks dramatically down when dealing with really big lists whose number is say 
10000 or more. Furthermore, people who know me won't be surprised that I also disagree with 
your abundant use of auxiliary functions as well as indexed variables. This may well be a matter 
of taste though. 

Anyway, would you believe that it is possible to write a program of this sort that doesn't contain 
a single indexed variable? Well, it is and here comes the proof! (Note that I assumed the matrix 
A to be sorted. If this is not the case, then use convert(sort(A)) as in the example below.) 

m_reduce_03 . dfw 

convert(A, A_ := [], i_, s_, t_) 

Loop 
i_ := 0 

s_ := FIRST(A) 
t_ := FIRST(s_) 

Loop 

#1: A :z REST(A) 

i_ :+ 1 
If A = [] 

RETURN REVERSE (AD JOIN (s_/i_, A_)) 

If FIRST(FIRST(A)) t t_ exit 
s_ :+ FIRST(A) 

A_ := ADJOIN (s_/i_, A_) 

#2: A := VECTOR([RANDOM(10) , RANDOM (10) ] , i, 1, 10000) 

(Simplify this line and delete the output before going on!) 

#3: convert (SORT (A)) 

Anyone who wants to try it out himself is advised to simply copy and paste the lines above into 
the input line. I did only a few tests, but it seems to work and its speed is awesome. Please tell 
me though, if you come across any errors. 

Cheers, 

Johann 

John Coleman 

Hello Johann, 

Your program is quite elegant and seems to be correct. It would be hard to break its speed record. I 
like your convention of using an underscore to indicate an identifier is a local variable rather than a 
genuine parameter. 

I am a bit puzzled by what you meant by performance breaking down dramatically with larger lists. 
With 10,000 element examples our programs seem to run about the same (slightly over 1 second in 
both cases on my machine in some runs, sometimes closer to 1.5 seconds with mine, depending on 
how close to sorted my ultimate output is). I tried A:= vector([random(10),random(10)],i,1 ,100000) and 
got the following results: 

convert (sort (A)) 21.1 seconds 

reduce_matrix(sort(A)) 33.0 seconds 
reduce_matrix(A) 60.2 seconds 

mine does seem 50% slower on large examples (in the sorted case to make them comparable). Is that 
what you meant by a dramatic breakdown - or did you observe even worse behavior? Just curious. 

Regards, 

-John Coleman 
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Alessandro Perotti 


per otti0 SCIENCE . UNITN . IT 


Dear Derivers, 

here is another efficient function which make the conversion asked by Marcelo Fernandez: 


m reduce 04.dfw 


Convert(A, a_) := 

If DIM(A) > 0 

#1: ADDOIN( [Ajljl, AVERAGE(a_ := SELECT(IDENTICAL?(u jl, Aild), u, A)’i2)L 

[] 


Convert(REST(A, DIM(a_)))) 

#2: Convert(SORT (A)) 

On big sample data (>50000) it seems (slightly) more efficient than other solutions. 

Best regards, 

Alessandro Perotti 

Johann Wiesenbauer 

Hello John, 

Thank you very much for your nice comments on my program. As to the "performance breaking 
down dramatically with larger lists", you are right, I should have been more precise. Below is a 
list of run times for various lengths k of the lists, where the lists I used are always of the form 

A := VECTOR([RANDOM(10) , RANDOM (10) ] , i, 1, k) 

The second and the third columns are the run times when approximating reduce_matrix(A) and 
convert(SORT(A)), respectively. 


k 

your version 

my version 

k 

your version 

my version 

10,000 

1.46 s 

1.17 s 

60,000 

108.50 s 

11.7 s 

20,000 

9.18 s 

2.88 s 

70,000 

171.5 s 

15.6 s 

30,000 

25.50 s 

4.52 s 

80,000 

223.7 s 

24.9 s 

40,000 

47.30 s 

6.55 s 

90,000 

284.1 s 

37.7 s 

50,000 

75.60 s 

8.48 s 

100,000 

300.2 s 

51.0 s 


Hence, I should have rather said that the "first signs" of a breakdown of your program become 
visible near 10,000, but the actual breakdown is only later. Interestingly enough the relative 
performance difference between our programs is biggest near 70,000 and becomes smaller af- 
ter that for some unknown reason. (Partly this may be due to the additional sorting of the list, 
which takes 7.62s for k= 100,000.) 

At any rate, I think that this has been a very fruitful discussion, in particular as regards the 
"power of indexfree programming". Thanks for your contributions. 


Cheers, 

Johann 
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Johann Wiesenbauer 


Hello Alessandro, 


Well, your program won't win a beauty contest, but it is admittedly the fastest one that has 
been been submitted so far as to this problem. 

Needless to say that performance is the top priority for me, hence it's my turn again. 


What about the following program that is at least as ugly, but seems to be again faster, e.g. 
4.21s for the case k=100,000 mentioned in my previous email! (Note that presorting A is no 
longer needed as well!) 


m_reduce_05 . dfw 

conv(A) := 

Prog 

A := REVERSE(TERMS(I(x a FIRST(v) • (FIRST (REST( v)) *y + z), v, A))) 

#1: If 3(FIRST(A), x) = 0 

A :z ADDOIN(FIRST (A) + FIRST(REST(A)) , REST(REST(A))) 

VECTOR ( [ x • 0( LN (v) , x), SUBST(v, [x, y, z], [1, 1, 0])/SUBST(v, [x, y, z], [1, 


0, 1])], v, A) 


|conv(A)| 


Cheers, 

Johann 


Alessandro Perotti 

Hello Johann (and all Derivers), 

now another throw by me. It's a small improvement (less than 30" on my PC for a 1 million 
unsorted table obtained approximating the line 

PROG (A : =VECTOR ( [RANDOM (10) , RANDOM (10) ],i, 1,1000000), []) 

m_reduce_0 6 . dfw 

#1: PROG (A := VECTOR ([RANDOM (10), RANDOM(IO) ] , i, 1, 1000000), []) 

# 2 : [] 

Convert(A, a_, b_, c_, t) := 

Prog 

a_ := KAPPENDCx 11 , v|l)-(v|2 + t) , v, A) 

#3: b_ := 3(a_, t) 

a_ := a_ - b_-t 
c_ := VARIAESLES(b_) 

[MAP_LI ST (DELETE (v, 1), v, c_), MAP_LIST(vilM2, v, [GRAD(a_, c_), GRAD(b_, O] 1 )] 1 

#4 : MWHA1 

Best regards, 

Alessandro 
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Johann Wiesenbauer 

Hello Alessandro, 


Well, I see that you are exploting basically the same idea, but using sort of indexed variables 
instead of powers. In fact, this works fine as to our current "test suite" for which it seems to 
be "tailormade". 


On the other hand, you have the serious restriction now that the first column of the matrix 
contains integers only as opposed to the other programs so far, which allowed even algebraic and 
transcedental numbers here. Hence, it's up to Marcelo now to decide whether he can live with 
this restriction or not. 


Cheers, 

Johann 


Francisco Marcelo Fernandez. Argentina 

Dear derivers: 

First of all, I want to thank you all for your interest in the problem of reducing a matrix with 
repeated entries. 

I had a practical application in mind and was having some difficulties with my awful pro- 
grams. When all the beautiful and fast proposals arrived, I had successfully prepared the experimental 
data for a further application of my (also poor) multi linear regression program using the first version 
of John Coleman's program for reducing the matrix. In my case the elements of the left colu mn are real 
numbers (which Derive dutifully converts into rational ones) and are commonly sorted. The elements 
of the right column are also real numbers. I apologize that the left column of my example contained 
only integers. Efficiency is not an issue, I am patient enough to wait for a few seconds to get the re- 
duced matrix. Besides, commonly the dimension of the experimental-data matrix is about 100-300 but 
my colleagues may eventually work harder and produce greater datasets. 

I want to size this opportunity to say that Derive New list is a remarkable help for those who 
are not good programmers and want to use Derive for practical purposes. 

Would you think I am a bore if I ask for a Derive program that makes linear regression with 
nonlinear parameters?. Mine is rather moody and decides not to converge more often than I like. 


Some days ago I met Johann at the University at the occasion of the traditional meeting of the Insti- 
tute. He regretted that he had not had time enough to compare calculation times of all presented ver- 
sions of this matrix reducing routine. Johann, this is my Christmas present for you. I compared all six 
versions reducing the same matrix consisting of firstly 50,000 and then 100,000 pairs of data. This is 
the result: 


Program 

m_reduce_01 
m_reduce_02 
m_reduce_03 
m_reduce_04 
m reduce 05 


50,000 pairs 
memory exhausted 
18.4 s 
9.81 s 
9.42 s 
2.52 s 
2.14 s 


100,000 pairs 
memory exhausted 
63.7 s 
26.1 s 
19.9 s 
5.17 s 
4.80 s 


m reduce 06 




It needs only little work with a graphics program to rectificate the pictures, Josef 






